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ABSTRACT 


A method was developed for forecasting swell using a 
Spectral wave-generation model based on a multiple point- 
mre Concept Ormsowelle=origin, The multiple point-source 
Goncept considers that the peak-energy swell emanating from 
Piovene cycloOnicestorm can be considered to have been pro- 
duced at one or more space-time point sources in the storm 
Baeeene impulsive antroduction of energy into the sea. The 


method was tested on five North Pacific storms generating 


/ 


swell recorded at Monterey, California. Predicted swell (Cgnce% 


ieehts which were made tometwo storms, were significantly 


lower than the observed heights. The time of occurrence of 


m@emprecdictea pedk height agreed with that observed for the 


—— = 


swell from one storm, but differed by about ten hours tO 





the other. Predictiens of the-dominant swell period were 


mecurate to-withipaapout one second over iliqueotiewer © lc 


i\ = 


Gov sca Vedeperiods for all storms. 
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le CBJECM EOS Thin oluU DY 


ine Ob jeetivewot thiswstidywas to develop a method for 
forecasting swell at a single observation site based on a 
Meta ple point-source concept of swell origin and using a 
spectral wave generation model. 

The multiple point-source concept of swell origin was 
mmeposed to this investigator by Professor Warren C. 
Thompson of the Department of Oceanography, Naval Postgrad- 
Meee oChOOlL. The concept involves the assumption that the 
peak-energy swell emanating from a moving cyclonic storm can 
Hemeonsidered to have been produced at one or more space- 
PiiempOolntc Sources in the storm by the impulsive introduc- 
tion of energy into the sea. 

The swell generation model proposed is a spectral model | 


utilizing the Pierson-Moskowitz (1964) spectrum and assumes 





that the spectrum of peak-energy waves present at a point- 
source can be computed from the surface wind speed at that 
point and the speed of the storm. 

ies products Of ties methodsare a forecast of the domi- 
nant swell] period and significant swell height with time at 
mime forecast station. 

The method was tested on five North Pacific storms and 
the forecast products were compared with the observed swell 


mecorded at Monterey, Calafornia. 





Hl. DEV EReR ENT ROS Tie SMe rnop 


A. BASIC APPROACH 

In a notable study of swell propagated over long distances 
aceerecorded offshore from San Clemente Island, California, 
Munk, et al. (1963) observed slanting ridges in the energy 
density topography drawn on a plot of frequency versus time. 
Each of these ridge lines described the frequency-time dis- 
PPrpution of the maximum energy in a swell train arriving at 
the recorder site from a given storm. It was shown that a 
trage line can be accounted for in terms of classical wave 
theory so long as the energy spectrum is assumed to have been 


Pee moitedaat~a point. source anspacceand timc. The slope and 





zero-frequency intercept of a particular slanting ridge line 
woiomenen used to compute the effective origin time and dis- 
momes trom the station to the source oi the peak energy 
swell. 

mice swe be pred ection modelwdeveloped an the present aeray 
effectively employs the reverse ot the above procedure. A 
"point-source" is identified in the wind field of a moving 
Syeronic Storm on each of a succession Of six-hourly weather 
maps covering the life of the storm. Each point-source 
represents a potential point of origin of the dominant wave 
energy emerging from the storm and is considered to be lo- 


cated at that point where the peak winds are generating waves 


an 


Shoewly “tOWdrd Stine (tot aiiieTatnomerOor Which the predierion 








is to be made. A wave spectrum is computed at each point- 


mource. the elardGteristiles, On CAChesSpectrum are determined 
from the surface wind speed at that point and the speed of 
mae storm toward the station. Each period component of the 


spectrum is propagated - toward the Station at a group ie Og 








icy which Seed Guin c tion Ofe1 ts | period. The dispersive ar- 


J 
ae on ee 


me a —— — — ee rca a 


rival at nee station of the waves Domenie CLs from i a single 


a tom ——— eee 


Pernt source -1s—nepresented bya plot of ner Ten Vier si) Ss time 








poet rival. 
The energy in each spectral component of the swell train 


represented by a particular ae Pane ebbing 1.6 earl cul ae 





by modif) modifying the energy spectrum at the. DO ly tes Onmecu OnmeriG 


effects of angular spreadi nidiehieeda —stelecitenGyaat tenuate 
SO HE ane FD tncfnat 


By graphing the energy-density values associated eh fore 
yar uous period components as a function of time of arrival 
Meeete Stat~on, anuvenergy density-time plot is obtained. 


Mitcmplot 15 termed a "propagated spectrum,. the envelope 





of the family of propagated spectra obtained from the series 
Papornt-sources during the lixe’ of a storm represents 
miiesenergy density of the spectral peaks at each arrival 
time at the observation station. 

tiemeeeotal ener oy ean othe swe ll vateanyactven arriValsstame 
is estimated by computing the area under the plot of fre- 
quency versus energy density values obtained from a time- 
mm or the propagated Spectra. The predicted (swell belents 
Pc SUC CEO DeCmS bails ML be lated LO thesic Goance 


Bmeregy GCalculatvons. The time distributvon of dominant 





erell periods 1S e0beaimed from the cime, distribution of (tle 
mrequency Components Containing the maximum energy density 
i the propagated spectra. ~ 

ie the £ollowmnc SccmmonsetNcscmpmacedures wil] be 


described and illustrated using synoptic weather events. 


eee GENERATION OF SEA SPECTRA 
ime voclection Of soca Level Pressure Charts 

The weather veiarts Used 1m this Study were six- 
hourly sea-level pressure analyses of the North Pacific 
Seecan produced by the Fleet Numerical Weather Central (FNWC), 
Momeerey. An advantage of using these computer-produced 
maivses 1S the objective manner in which the isobaric fields 
are produced. The fields analyzed covered the entire area of 
the North Pacific under investigation. 

The weather charts were chosen on the basis of 
frequency-time analyses of five selected swell trains 
mecorded at Monterey. Each of these swell trains had been 
analyzed for height and period of the dominant waves. The 
eorm producing CaGieeet der doe CCIE ie ¢ GmOnutine wea Emeds 
maps from the effective spesagn JENS BHO. ll SieehNelS a eneal eval 
from the observed frequency-time distribution in the swell 
mimene Manner described by Munk; et al. (1963). All were 
North Pacific storms. The sea-level pressure charts cover- 


mie the important part of the history of each storm are 


Waisted in Table 1, along with the corresponding swell trains 


recorded at Monterey. 





Table I: STORM AND SWELL DATA 


A. North PacGittre Ste:ms otudied 


storm 
Designation Sora CCme Me soumemciduauSm ma hy. Zed e( GME) 
1 W200/15 Fep to 200/17 Reb 67 
2 Le Ii Vreie sere) Le Al stele (or 
3 1200/12 Nov to 1200/15 Nov 67 
4 0000/16 Nov to 0600/19 Nov 67 
5 1200/12 Mar to 0600/15 Mar 69 


B. Associated Swell Trains Recorded at Monterey 


swell Train 
Designation Observed Swell Occurrence at Monterey (GMT) 


iE 1200/20 Feb to 1600/23 Feb 67 
y 0800/23 Feb to 0800 25 Feb 67 
5 PA007 LyaNeveron 2200/7 ZU Noy ov 
4 0400/21 Nov to 0600/23 Nov 67 


S 2000/18 Mar to 0800/20 Mar 69 


10 





2. Determination of Surface Winds 
By Mé€asuring the spacing between isobars on sea-level 


Bressure charts the average geostrophic wind speed was cal- 


culated from: ve K AP 
9 22NP sing An 
where ee geostrophic wind speed in knots 


Ree= pressure ditterential in mvliibars 
Oo Sobanespacinewan degrees latitude 

fi SelevcalwinclsS shih Wewveee 

P= (eisai tay Oteeea ad. oer le Cee oles om/cm>) 


N= Uiitseeonvers@ton- cacror (1.7725ex% i9°° 


) 


Q = earth's angular speed of rotation (7.29 x 10° 
radians/sec) 


5 


In this study, because of the relatively small scale 
Sietme weather maps used (1:60,000,000), geostrophic winds 
were calculated using a 12-millibar interval (3 isobar spac- 
ings). Geostrophic winds were converted to surface winds 
Peseta Surtace to geostrophic wind-speed ratio of 0.8. 
FNWC currently uses an average PatronoL Uys based on 
Carstensen's (1967) report showing the variation in the ratio 
with latitude (Figure 1). In the present investigation, the 
five storms studied were located between latitudes 31 and 
m> degrees North; accordingly, a ratio of 0.8 was considered 
mombe a more representative average for these latitudes. 
The surface wind direction was assumed to be at a cross- 
movdr angle or Usedecrees; = tnis agrees with current —Fllect 


Numerical practice. 


ik 





FNWC 
average peed 


North Pacific 
average 


(deg) 


Latitude 





Figure 1: Surface Wind-Geostrophic Wind Ratio (Vs/Vg) as a 
Function of Latitude (after Carstensen, 1967) 
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Se Location cl athe mies t -oeunee 

From consideration of the isobar curvature and 
surface wind-speed profile characteristic of a cyclonic 
storm, it is assumed that waves of maximum energy propagat- 
mime toward Monterey must be generated at or close to the 
point where the maximum surface winds are directed toward 
hemcerey . 

li Order stemudcitueny a enat location in the cyclonic 
wind field on a selected weather map, an overlay of great- 
Mmirete arcs from Monterey was prepared (Figure 2). By 
moat ine the overlay fifteen degrees clockwise about the 
aporm center the locus of points in the storm where the 
meeemercircie arcs are tangent to the isobars can be deter- 
mined; this locus represents the storm radius along which 
are found those surface winds generating waves directly 
toward Monterey. The point of maximum wind speed along this 
mcs 1S identified by measuring Successive 3 isobar 
(12-millibar) spacings and locating the midpoint ot tne 
Benoengest pressure gradient. Figure 2 illustrates this 


meenmrque Lor a selected surface pressure chart. 


AL SVeetr Teac lon of the Sea spectrum at the Point-Source 


a. Basic Considerations 
Moskowitz (1963) analyzed wave records and 
corresponding wind measurements taken by ocean weather ships 
of the United Kingdom in the North Atlantic Ocean and 
ebtained a family of spectra representing fully-arisen sea 


MoOnadlt Vols tar Certaan wind speeds between Z0 and 40 knots. 


14 





Pierson and Moskowitz (1964) later used the data to test the 
Similarity theory of Kitaigorodskii who proposed that all 
fully developed seas should have the same non-dimensional 
Spectral form. They showed that within the limit of accu- 
racy of the reported wind speeds, the data of Moskowitz 
@eteirmed the theory that a single non-dimensional power 
meer um Might be used fo Fepresent the dimensional spectra 
for fully arisen seas generated by any wind speed. The non- 
faemensional spectrum reduced from these data was determined 
Mommave 1tS non-dimensional peak value of 2.75 x 1 = at a 
non-dimensional frequency of 0.140 using the following non- 


fimenslLOonalizing conversions: 


v= f V. /g | (1) 
~--- 5 : : } 
S(5} = Elf) 9// Vv, (2) 
where f = dimensional frequency in hertz 


f = dimensionless frequency 


E(f)}= dimensional energy density at frequency f 
_§(¢) = dimensionless energy density at frequency § 
V,; = surface wind speed in m/sec 


Uee- ACcelehkatlonnon gravity (9. SUSS m/sec¢) 
PrCuschncicc CLO me MOntCameniat the general 
analytical formula for a one-dimensional gravity wave spec- 


erum can be written: 


Sf ake tee ot) (3) 


Ws 





where S(f} = spectral energy density for frequency f , and 


where it can be shown that: 


co 
T 
3/3 
eo 
s 


© mn 
K,= f, S(f,Je 
{oe See ici es aeOremtom= De GUrd lmpedi 


[iit eee Weomc mb retseanelder (1963) proposed»d 
Spectral form for wand-pgenerated seas based on Equation 3 in 
which m = § and n = 4. Pierson and Moskowitz concluded that 
fmommeince purpose of Lorecasting the properties of the larger 
Mawes itl a Wwind-generated sea, this spectral form is a good 
memnesentation. 

Dee GCenemrtroOne On bile ey opcet ra 

In this investigation, Bretschneider's analytical 
formula (Equation 3) was used to generate Bretschneider's 
(1963) spectral form approximation Oa tile tems On -MoOSs koh 
non-dimensional spectrum over the dimensionless frequency 
manee 0.10 to aot eusingeans  BMeVoZOcenputer,  Fipure 3 is 
pepplot of this spectrum. 

The dimensionless spectrum was used to obtain a 
dimensional power spectrum at the point-source on each 
fetiver Ghart wsing the wind speed value calculated for the 
point-source and the conversion factors of Pierson and 


Moskowitz (Equations 1 and 2). 


AG 





-5 


x 10 


Energy Density 


Dimensionless 





Bieure $: 


Dimensionless Frequency 


Dimes Tomress Ove specu num 
(after Bretschneider, 1963) 


de 





The spectrum calculated using this procedure 
mepresents a2 fully arisen sea veonditiom. Since the sea at 
the point-source on most weather maps was clearly not fully 
moesen, the tulily arisemespectrum was truncated at an appro- 
beaate cutolLft frequency, as proposed by Neumann (1953), to 
=eproximate the spectrum of the non-fully developed sea con- 
sidered to be present. The basis for the procedure devised 
follows arguments presented by Braunstein (1970). 


Braunstein aepied that fully arisen seas should 


ea a a 
4 





Demgenerated when the velocity of the fetch equals ithe group 


merecity of the peak Swell being generated, since such a 


= 


ze = ———— 


meencrdence produces an ettéectavely one duration. sin the 


gare enamine 2mm ~~. 


a ee 


peescent investapation a group velocity-fetch velocity 








een 


—_— 


Bre racnce was used to crteqiaa -duration imat,onethe fully 


eee penne et TE AIT 


arisen enerey spectrum. If one travels with the dominant 
waves at nee Ceelpwe LOGIE seUlcmecne Gy apt mItT tO ticused 
meeolven wand Speed 2s a function of the duration over 
which the waves remain under the influence of the generating 
wind field. In the case where the waves travel in the direc- 
tion of storm movement, as is the circumstance in this study, 
fae longest effective wind duration is achieved for that 
frequency component whose group velocity equals the speed 

of movement of the generating source. Lower frequency com- 
momentS, 2f generated, move faster than the generating source 
ano move ahead of the area of peak winds before they acquire 
full energy. Higher frequency components require a rela- 
tively short wind duration for achieving saturation and 


Gemain in that state. 
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In this study, Neumann's (1953) model of a 
iar p iS pectralmecltOrimmas ioclned amd the Spcctrum gener- 
meed by the peak winds at Che point-source on a given 
weather map was truncated at that frequency having a group 
molecity caval to the Speed of the point-source averaged 
ever the six-hour period just preceding. 

Tiss netnodshasssome Obvious deficiencies. For 
feuy large storm speeds, 1t tends to produce spectra which 
are more developed than would actually be the case, and for 
very low-speed storms (including the stationary storm case) 
fiicmspectra proauced are less developed than would be 
@eepected. As will be shown, however, for forecasting the 
characteristics of the dominant swell produced in a moving 
Seronle Storm, the method produces results which are 
reasonably consistent with observations. 

Dice aortei Oumspccrl yim tnewdtmattOn assoc? 
meeaswith the peneration of the maximum seas in a moving 
eyerornic storm of given Size, intensity, and speed is not 


simple and merits: additional research. 


Sree ONELL PROPAGATION 


uo per sion or Period Components 
Thewspecthal components of the predicted seas are 
Bich permitted Lo propagate to Monterey. The arrival time 
at Monterey from a given point-source is computed for each 


component using the dispersion relation: 


ie ef et 4 dg 


Ie, 





arrival’ Cine 01 speri0d component 1 


I 


where t 

f = origin time at the point-source (a weather-map 
time ) 

de moredteci Clewdistance trom Ene POInt-source to 


Monterey 


T = period in seconds 


Figures 4 through 8 are plots of swell period versus 
mame O£ arrival at Monterey for each of the five storms 
meuaied. For each storm there is a set of curves, each of 
Demtcm represents the distribution with time of the compon- 
ents of the wave spectrum arriving at Monterey from a given 
point-source. Each curve is labeled Witha date-time group 
indicating the map time of the particular sea-level pressure 
Sime CONtCaining the point-source from which the curve was 
womeed,  iie tength Of “each Curve 1S poverned by the char- 
feucristics of the sea spectrum calculated at the originat- 
ing point-source. The high frequency end (that component 
wariving latest in time) is determined by introducing a 
ioxed Minimum energy density criterion into the dispersion 
program; frequencies inthe sea spectrum having energy den- 
Sity values less than 1.0 m?-sec Phemnomplottea, i Ihewlow 
frequency end of each curve represents the frequency of 
epectral cutotft discussed in the previous section. 

Zo Prom aea honor oe Cen d lame toy. 
a. General Considerations 
DiadcibetoOi tOrunewerFeCt OL dispersion, the 


ener ey In thie wave specunum Scienated at each point-source 
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ms diminished by the effects ox angular spreading and high 
frequency attenuation during propagation. the angular 
spreading function is the frequency-independent Neumann 
function (Figure 9). An attenuation function is used to 
account for the non-linear processes of wave-wave interac- 
tion and wave breaking; it results in dissipation of energy 
mainly in the high-frequency end of the spectrum with little 
noticeable effect on the low frequencies. 

b. Angular Spreading 

Although on each weather map the maximum energy 
waves being produced are considered to have originated at a 
point-source, the generating area has a finite width by 
virtue of the fact that the maximum surface wind speed com- 
puted at the point-source is obtained by averaging the wind 
meena 12-millibar isobar interval. The effective width of 
Mes ceierating area 1s therefore taken to be equal to this 
12-millibar spacing. 

An angular spreading factor is calculated using 
the fetch width, the great circle distance from the point- 
source to Monterey, and the trigonometric relationships 
outlined in Pierson, Neumann, and fenes (1955 es. ius: 
trated in Figure 9. Because the surface wind at the point- 
source is directed toward Monterey, the effective fetch 
width as viewed from Monterey has a fixed angular relation- 
ship to the great circle trajectory from the point-source 
to Monterey. This circumstance simplifies the geometry of 


angular spreading and, along with the adoption of a fetch 


26 





waath defined by the@pressure exadiemt, makes it possible 
to calculate angular spreading factors as a function of 
@rstance and wind velocity. 

For practical use, a table of angular spreading 
factors was computed for wind speeds of 10 to 70 knots, at 
merements of great circle distance of 60 nautical miles 
(20r distances of 600 to 5000 NM), and at latitude incre- 
ments=of-one -deerce ¢from 30 eto 70 degrees North). 

In the present study the swell are propagated 
eemeracally along great circle paths; accordingly, the PNJ 
method, which assumes straight-line propagation, was modi- 
ted using spherical trignometric relations. <A comparison 
of the calculations using both methods indicated that the 
mermene Increase in accuracy obtained using Spherical 
eraconometry was not sufficient for the propagation ranges 
considered (2500 to 5000 nautical miles) to warrant the more 
Boupli1cated spherical computations. 

CP hon ee ney mrt cilia tron 

Although the subject of wave-wave interaction 
miamether non-linear processes has received considerable 
treatment in the recent literature (e.g., Hasselman, (1963), 
Phillips, (1963), and others), no analytical representation 
ireepeen derived that 1s suitable for application to long- 
range swell prediction. However, Snodgrass, et al. (1966) 
Maccoented empirical data on the attenuation of selected ifre- 
mmencles Observed in their study of long-range propagation 


me OGeal swellwin the Pacatic. in general, they found 


a7 


attenuation to be large within the limits of the wind area 
Seethe generating storm and negligible beyond the storm. 

mr swell that had travelled a considerable distance from the 
eeorm, they also found that below a frequency of 0.06 hertz, 
attenuation was too small to be measured and above 0.08 
hertz, the wave energy was masked by the background level. 
The empirical attenuation data given by Snodgrass, et al. 
omer lorvarithimic attenuation coefficients reported in units 
of decibels per latitude degree of propagation distance. 
Merethne range of fréquencies 0.06 to 0.08 Hz, these data fit 
an energy attenuation function of the form: 


—-2ax 
e 


iT 


where a@= modulus of amplitude decay in deg 1 =O .d dts 


Ke.) 
iH 


logarithmic attenuation coefficient in dB/degree 


ba] 
i 


Propagation distance in degrees 


In this study, the empirical data of Snodgrass, 
eer. were used, and were extrapolated to provide an e€sti- 
mate of the most likely attenuation at higher frequencies. 
ne extrapolation of the aLeenuatvomedaata LSmmiusStraved an 
Figure 10, a semi-log plot of the logarithmic attenuation 
coefficient in decibels Pel elatnbpiecmocOrCemuCr SUS Lrequcheg, 
Gime empirical data were extrapolated to the frequency at 
which the logarithmic attenuation function reached a value 
of 1.0 dB/degree; higher frequencies were assumed to be 
equally affected by attenuation. Frequencies below 0.06 
hertz were similarly assumed to be attenuated by the same 


amount as the 0.06 hertz component. 
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Figure 10: High Frequency Attenuation Function 


50 





wes 


8 Storm Limit 





Monterey 
- © 


Figure 11: Storm Limit for Surface Pressure Chart 
0000Z/13 Nov 1967 


31 





Following the oObSeryvapronvot Ssnede ris > ec cmc. 
that attenuation beyond the area encompassing the storm fetch 
was negligible, the energy attenuation function of Snodgrass, 
et al. (1966) was applied over the propagation distance from 
mie point-source to the outer limit of the storm. The storm 
limit was taken to be the boundary beyond which the isobars 
momronger follow the circular pattern of the storm and thus 
Saneno~longer=bewassocrated with the counter-clockwise flow 
around the low-pressure system. Figure 11 shows the location 
Seethe storm limit on a selected sea-level pressure Chart. 

d. Propagation of Spectral Energy 

inewenere, dissipating effects of angular spree. 
ing and attenuation during propagation of the swell from the 
femme SOUrCe On Cach successive weather chart to Monterey is 
@ecoumted for by applying the spreading and attenuation fac- 
Mere, to the energy-density values of selected periods cover- 
Mieerhe entire spectrum. For a given frequency component, 


ie clergy density propagation loss is computed by: 
| 2ax 
a (4,) = [ale | Le 


energy density in deep water at Monterey 


i. a) 


EG = energy density at point-source 
Fe = angular spreading factor 
a2cx ; 
e = Elie EINE a ooueTbigle re Mane 


By propagating each of these dissipated energy-density 
maolues to Monterey at the appropriate group velocity, the 


distribution of spectral energy with arrival time is obtained 


OZ 





mer C€ach point-source. jis enero Scdeusily-cime curve wis 
meserred to here as a “propagased speerrum = Plots wot tic 
propagated spectra computed for deep water off Monterey for 
the five storms studied are shown in Figures 12 through 16. 

In each figure the uppermost envelope that can 
be drawn to the family of spectra shown is considered to 
represent the total propagated spectrum for the actual 
swell train that would be expected off Monterey. It may be 
@eeneain the case of Storm 1 (Figure 12), that the spectrum 
generated at the time of the weather chart of 0000Z/16 
February 1967 had sufficient energy in all frequencies to 
aoe te the other spectra arriving off Monterey. For the 
other storms, it is seen that the swell at Monterey is 
omamated at different times by waves arriving from differ- 
pumeelocations along the storm path. 

The frequency associated with a particular 
energy density in a selected propagated spectrum cannot be 
read from Figures 12 through 16, but may be determined from 
the corresponding period-time curve (Figures 4 through 8). 
fwereeoectra shown in Figures 1Z through 16 are greatly 
mene d iil CMS Tin TINS lea Ie eerie Ss Sieve ce) Telave 
How trequencies compared to the spectra present at the 
Beurce-points. The truncated propagated spectra that may 


be seen in the figures are non-fully arisen spectra. 
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Dee PREDICTION OF SWELL ChARACTIERI ST ICS 

Ie Jae Welecoe Jeno sore Ieee seis aber 

Each curve in Figures 4 through 8 potentially 
meapResecnts the period-time Variation of the dominant energy 
swell arriving at Monterey from one of the five storms. The 
Pamirecular curve or portion thereof containing the dominant 
wiemey at ad Piven moment Can be determined by referring to 
femmes 1Z throtigh 16. The upper envelope of the family of 
Meeepacated spectra Shown in each figure represents the time 
mestraibution of the dominant energy waves in the combined 
ier arriving from all point-sources. The date-time group 
Pimirlryine Cac Spectrum rerers to the map time of the 
pomme-SOurce from which the spectrum was derived, and can 
Pemvscasco relate the spectrum to a particular period-time 
@mewe (in Figures 4 through 8). 

Mice PeErLoOnsceOnstNesemmerilod time Gurves containing 
wiemaominant energy were 1déntitied in this way and appear 
Mteures 17 through 21 as heavy lines. The dots plotted 
Mmeten figure represent observed period data and are dis- 
aesed In a Later section. 

. Rice Gergen onOulenotemrnadt ther ticavy lines 1n 
Peines 1/ through Zl represent the predicted distribution 
with time of the dominant swell periods in deep water at 
me sSclson Site. | ihe eftects on Shoal-water modifications 


to the arriving swell are discussed later. 
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2. Deep Water Wei ght Prediction 


a. Estimation of Arriving Sie Energy 

A moving storm may be viewed as the source of 
mimincinite number of energy spectra Penen aed lobe ipsa 1c Ve 
meor S lifetime. the family of enerpy-time curves for 
Pomeerey shown in Figures 12 through 16, then, are derived 
from a finite number of selected spectra; specifically, one 
ppeetrum associated with each weather map at the point- 
memrce identified. By making a time-cut through one of 
these composites, the energy densities of the spectral 
Gomponents arriving at Monterey at that particular time can 
be read from the propagated spectra. The period of the 
waves associated with each of these energy-densities can 
be read from the corresponding period-time curve (Figures 4 
tnrough 8) by making an identical time-cut. Thus, for any 
Sewmeeted time, a plot can be made of energy density versus 
frequency. A continuous spectrum drawn through the points 
Of such a plot represents the predicted energy spectrum of 
meemeswell in deep.water at Monterey at that time. Figure 22 
metstrates these predicted spectra for a series of arrival 
times for Swell Train 3. The area under each spectrum 1S 
feeportional to the predicted energy in the arriving swell. 

bes Calculation of Swell Heights from Energy Estinaees 

ihe total enercy ino asea tomeeund by atc cearee 
ing the energy spectrum over the full range of frequencies 
Peon WW tO co, This energy value, 1 the approximate opece 
tral form is known, can also be related to the statistical 


distribution of wave heights in the sea. Longuet-Higgins 


4S 


(1952) has shown that for narrow band spectra (i.e., swell) 
the significant wave height, the average of the highest 


poerd of the waves in a wave record, 1s piven by: 


= 2. E 
Hy, 2.831 (4) 
where sh significant wave height 
3 
E = total energy in the spectrum 


If the predicted swell spectra are sufficiently narrow 
Banded in the sense implied by Longuet-Higgins (1952), 
01S See Great relationship provides a means for computing 
the swell height distribution with time from the predicted 
meme y-time distribution. 
chet pdetuoneaned Shoaling Modrtications 
an SIOCichCh wc] OpmemE Gedue & ron 

The wave sensor used to obtain the observed 
eels height at Monterey for the five storms studied was 
ieeated in shoal water; accordingly, the effects of refrac- 
tron and shoaling had to be considered in order to compare 
observed heights with heights derived from the computed 
swell energy. 

Simcee retraction effects are dependent on the 
direction of swell arrival, and directional information was 
Momerecorded, two procedural alternatives presented t#hem- 
eelves: the arrival direction of the observed swell could 
be estimated, and the observed heights adjusted to remove 
micwcitieets OL Tetractlon anu sShodling: on tne predicted 


heights could be modified for these effects and compared 
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Wath the shoal water data as observed ab Cteeccn or wo ace 
The latter course of action was chosen because an arrival 
mzimuth for the swell werain propacated trom cach peine- 
meurce 1S a Lunction of eine great circle path between the 
point-source and Monterey. 

To. Tllus Clave sinew: Cecadume Used mre rca elma. 
meemrecioht distribution in the swell at ytne wave-cape sire. 
Seorm 3 will be used as an example (Figure 22). The 
individual energy-density values determined from the deep 
feeeer Spectral time-cuts, indicated by the dots in the 
figure, are modified for the effects of refraction and 


iealing using the relation: 


fatal] = fastid] Ke K. 1 (Ss) 


where asst] shoal-water energy density associated 
Wanteie 1: cequeney at 
[asc] = deep-water energy density 
K. = refraction coefficient : 
K,; = shoaling coefficient (water depth 30 
| ei) 
Refraction coefficients for the sensor site are shown in 
Peeure 23; We data were compiled for the sensor site from 
a number of manually drawn refraction diagrams prepared as 
mmlaboratory exercise by fLormer Students at the Naval 
Rost egraduate School. 
TMEwPeENer ey Inmseme Swell arraving at the sensor 

Site at a selected time is then approximated using the modi- 
fred energy density values in the step-wise integration 


scheme: 


4/7 





n-—j 


5st) =} [Acted] 3 Pathe] fF (6) 


p= 


minere E(f) = total shoal-water energy in swell arriving at 
Game sit 
2 
lA,(s) = modified energy-density value associated with 
frequency f. 
n = number of energy-density values used in the 
approximation 
Af = frequency band between f, and § 


i+ 


Figure 24 illustrates the shoal water energy spectra for 
sell Train 3 at selected times. 

Predicted values of significant swell height in 
Shoal water may then be computed from these shoal-water 
Sreeroy computations using Equation 4. Predictions of signi- 
ficant wave height for the swell arriving at the wave gage 
ere trom Storm 3 are shown in Figure 26. 

baolOdl Water Perrod Predietnon 

NSuaeGesultwOt LTetracrion andeshoaling.. tne 
energy density associated with a given period may be expected 
to differ at a shoal water site from the deep water value. 
Peeause these processes are frequency dependent, a shift of 
the iGadinesnt Chieti) ila) POC EErONmonesocrrToc- time Curve to 
mrotner aS the swell passes from deep to shoal water. 

Pe prcche Ol Ome mCi anmemmCrTOGS Ln Shoal 
water is obtained by plotting the periods of maximum energy 


density in the shoal water spectra as a function of arrival 
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fone. Figure 28 i2lustrates Ghevpredictcd (91) 0 aime Ween 
meniod-time curve tor swell lraam 3 asede temminedst rome 


Shoal water energy spectra shown in Figure 24. 
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Li tvoncGOMPARISON OFePREDI CUED SAND SOB SiR) Piel Eee 


A. OBSERVED WAVE DATA 

The observed wave data used in this study for comparison 
with the computed wave data were obtained from analysis of 
frrip-Chart recordings from two bottom-mounted pressure-type 
Peemocensors located 1n approximately 30 feet of water. Hie 
Strip-chart records from both sensors were manually analyzed 
iy sEroressor Warren C. Thompson. Wave periods were obtained 
using the wave-group method, which for distant swell has 
been shown to yield periods equivalent to the periods of 
maximum energy density obtained by spectral analysis 
(Thompson, 1973). The wave records exhibiting the swell 
Smucinating from Storms 1 through 4 were obtained from the 
Monterey sensor located at Del Monte Beach. The swell 
periods from Storm 5S were obtained from a sensor at Stinson 
meaen, California. The Stinson Beach sensor 1s located 
about 75 nautical. miles closer to the swell source than the 
Memeterey sensor and approximately on the same great circle 
i~amectOryesoince the predicted periods weregcomputed for 
Monterey, the arrival times of the observed periods 
recorded at Stinson Beach were adjusted to correspond to 
arrival times at the Del Monte Beach site. 

Significant wave heights at the Monterey gage were 
determined by adjusting the apparent heights recorded by 


the pressure gage for the effects of hydrodynamic damping 


Se 





due to water depth using Wiegel's Tables (Wiegel, 1954). 
The recorded swell was presumed to be narrow banded so that 
a single damping factor could be applied with good approxi- 
mation. The pressure-response factor of the wave gage was 
greater than 0.988 over the period range from 5 to 20 
seconds. The significant wave heights computed from the 


Pmessure gage data are reterred to as observed heights. 


Be OWELL HEIGHT 

Significant swell heights at the Monterey gage site 
were predicted for the swell arriving from two of the storms 
Studied (Storm 1 and Storm 3). Comparisons of predicted and 
G@bserved signficant heights are shown in Figures 25 and 26. 
It may be seen that in both sets of swell the observed 
mementS are ¢réater than the qimedicted heights; in Swell 
Memin | the observed heights exceed those predicted by a 
maetor of about 2.6, whereas in Swell Train 3 the factor is 
Titties lneeiemcase Of Swell Traam 1.8 the opserved peas 
Mement Occurred approximately ten hours earlier than jee= 
emetcd. Ihe peak height in Swell irain 3 occurred approxi- 
Hately at the time predicted. 

Wie dtirerence 11 Negcmitude Dbetweem eneuprearcted and 
Spserved heights most likely can be attributed to one or 
more of the procedures used to apply angular spreading, 
mecenuation, and shoaling and refraction. Ihe choice of “a 
m-millibar imterval tor demining ther rerenmwirdti1 15) langeds, 
arbitrary and may be too conservative. The way in which 


the storm limit is defined may have introduced excessive 
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Figure 28 





attenuation estimates. Poor refraction data also may have 
mmteroduced a Sieniticantscrnor wn the ca heuterloncwoceenc 


enoOal-water swell energy. 


C. DOMINANT SWELL PERIOD 

The distribution with time of the dominant swell period 
at the wave gage site also was computed for Swell Trains 1 
gma 5. the results are shown in Pigures 2/7 and 28 by open 
mecles; the observed periods are represented by dots =) ihe 
mmecaicted and observed periods at the sensor site are supen- 
moposed on the family of potential period-time curves com- 
muced for Monterey as shown in Figures 4 and 6. As stated 
mameier, the heavy portions of the perrved-time curves 
mepresent the predicted dominant periods in deep water at 
Momlterey . 

in Figure 27 at may be observed, from the predicted 
period distribution with time in deep water and at the wave 
meme Site, that the shoaling and refraction modifications 
Mme propagated Spectra resulted in a Shift of the doman- 
emt energy trom - DerTod tiles Cll ye hom Mlereiicit lam sec 
period bandeirom 18 to 2) seconds. ~Filpupe 72 seexhpmprtsea 
aouilar Shalit: however, in this case the predicted shift is 
more consistent with the trend of the observed periods. It 
Mey be concluded that thesdomimant spermodmean shiit due to 
Emoaling and relLraction, “bUG wnatathe soiree ou mtOt lance. 
mecovding ily, the predieteumperiod—-tinerd met i buttons in 
deep water are used for comparison with the observed peri- 


ods at the wave gage for Storms 2, 4, and 5. 
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Comparison of thewperiod= tines diSetl abut 1 ono see eit 
dicted dominant swell (heavy lines) with the mean trend of 
mac observed periods (dots) for each of the storms (Figures 
17 through 21) shows close agreement. For the five storms 
memerecd, the mean difference between the predicted and 
Meeerved periods does not exceed about one second. It may 
memeoncluded that the method produces accurate period 
ieerecasts. 

Pt Cerences mien wewicl ome ratelimne  Cimdmacte hist les 
among the five storms may be noted from examination of the 
ferecicted periods (heavy lines) and the observed periods 
Heocs) shown in Figures 17 through 21. In Storm ] (Figure 
Mme, the period-time curve propagated from a4 single point- 
source (on the weather chart of 0000Z/16 February 1967) 
somcains the dominant energy over the entire duration of the 
meererved swell at Monterey; this is veritied by the close 
memnespondence of the predicted deep-water period distribu- 
mace with the ebserved periods. In Storm 4 (Figure 20), 
Mm@iceaistribution of observed periods also indicates effec- 
Mee CCheration at a single point-source; however, predic- 
tion called for the longest period components of the swell 
mompe Senerdted at the point-source identified on the 
mmeceecding weather lap, and )tOudrrive sooner. it 1S evident 
feo an e€Xamination Of themplot ot propagated spectra 
@erived £Lrom thas stemm (Faorre ! 5 )5. that qercuminatiomn ol 
Wiech PELIOd= Line Guiness meeMi ami eett eco tne tlic een 


Enis Case particularly dependent upon the accuracy of the 


Oe 


modifications to the sea epee (num escaCchepoint soumec 
due to spectral cutofEs aneular sprmeadine, and attenuacien. 
By CONCTrast, TM Sess eee cer lo eet o rem Sac dinec. 
moe Shift of the predicted dominant periods as well ads the 
observed periods from one period-time curve to another such 
immet the dominant energy arriving at different times is 
meomagated from several different point-sources. In Storm Z 
(PassuresdS) ,waecdmabar shiftiin.the dominant periods across 
me, GiSpersion Curves may be noted. In this storm, the 
mervod-time curves are nearly coincident, which is due lto 
the fact that the speed of the point-source toward Monterey 
and the group velocity of the peak energy component in each 
sem ce spectrum were closely coincident. The wide scatter 
@i periods recorded after 0000Z2/25 February is attributed 


mona Wave train £rom another disturbance. 


eee REFINED PERIOD PREDICTION 

I OLGCCETEOg un ther 1ivestreoatcmtic sap aie tarsi ft on 
oe dominant energy from one pemod-time curve to another, 
which 1s particularly apparent in the case of Storm 3 
moroure J19)> the synoptic parameters derived fe One d Gla sli 
fatee pressure chart (i.e., surface cand Speed. Storm Speed, 
and distance from the point-source to Monterey) were 
miecarly iInterpotated dt two NoMneaiercnicness» lhe yperlod- 
eime distribution Of the dominant svelmeae Monterey Obtained 
from this refined data base is illustrated in Figure 29 by 
wae heavy Jines Ihe assocratred propacatea energy Speeund 


are illustrated in Figure 30. The reader may note that the 
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time scecquence of the dominant segments on the period: time 
eurves reveals the generation ynastory. ot thevarriy ine vors 
erain. In the period band from 18 to 25 seconds, the domin- 
mime energy was generated am the later part of the storm 
history (1000Z/13 November through 2000Z/13 November); the 
menmrtod band from 13 to 18 seconds was generated early in the 
storm history (12002Z/12 November through 1800Z/12 November). 
ine -cencralwtiwend Of thesdeminant periods in the pre- 
feed ca swell may be Seen tO appuoximate the trend of the 
@escerved periods; however, the correlation does not appear 
boemaye been significantly improved by refinement of the 


e-hourly synoptie data. 
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IV, CONCLUSIONS 


In spite of the poor height predictions obtained for 
the two swell trains at the Monterey wave gage site, the 
Tebeiple point-source methodwis surtrerently promising with 
tiereard to the time of arrival of the maximum swell height 
and the accuracy of the dominant periods obtained that 
additional research effort should be expended to test the 
iethnod further, and to tune the procedures involved so as 
mmoive better height) predictions. The latter may be 
accomplished by modifying some of the procedures within 
fealistic limits, particularly the procedures for angular 
Beereading and attenuation. W8ie predicted heights should 
be compared with observed wave data of better quality than 
were available in this study, preferably recorded in deep 
water so that the complications of shoaling and refraction 


can be eliminated. 
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